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Reconciling Planck results with low redshift astronomical measurements 
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We show that emerging tension between the direct astronomical measurements at low redshifts 
and cosmological parameters deduced from the Planck measurements of the CMB anisotropies can 
be alleviated if the dark matter consists of two fractions, stable part being dominant and a smaller 
unstable fraction. The latter constitutes ~ 10 per cent at the recombination epoch if decays by now. 


INTRODUCTION 

Measurements of the cosmic microwave background 
(CMB) fluctuations by WMAP pQ and Planck [2 Col¬ 
laborations opened a new era in high precision cosmol¬ 
ogy. These data are well described by the standard 
spatially-flat ACDM cosmology with a power low spec¬ 
trum of adiabatic scalar perturbations, and make a great 
step towards the precise determination of cosmological 
parameters, in particular of the Hubble constant Hq = 
100 h km s -1 Mpc -1 , the density parameters f lb and Uldm 
for the baryon and dark matter fractions, and therefore 
of the whole matter density Cl m = fl b + Uldm = 1 — Da- 
Recently the latest results of the Planck Collaboration 
have been published [3] based on the full mission Planck 
data. They are in excellent agreement with the 2013 data 
[2] but with improved precision. 

The Plank data imply a rather low value for the Hubble 
constant, which is in tension with the direct astronomical 
measurements of h. The Planck 2015 TT,TE,EE+lowP 
data, which we take in this Letter as a benchmark, de¬ 
termine the Hubble constant with 1 per cent precision, 
h = 0.6727 ± 0.0066 0. 

Direct astronomical measurements of the Hubble con¬ 
stant indicate larger values. The analysis j4j of the 
Hubble Space Telescope (HST) data based on over 600 
cepheids in host galaxies and 8 samples of SNe la yields 
h = 0.738 ± 0.024 including both statistical and system¬ 
atic errors. Independent analysis of Carnegie Hubble pro¬ 
gram [5] using the Spitzer Space Telescope data for cal¬ 
ibration purposes lead to h = 0.743 ± 0.021. Both these 
results are discordant with the Planck result at about 
2.5er level. Other astronomical estimates also typically 
imply high values of the Hubble constant. For exam¬ 
ple, the analysis of the gravitational lensing time de¬ 
lay measurements of the system RXJ1131-1231 implies 
h = 0.787 ±4.5 0. 

In addition to h, there is tension between other CMB 
derived observables and their direct low redshift measure¬ 
ments. The Planck results show a tension between the 
cosmological constraints on as and fl m from the CMB 
[3 El and from clusters as cosmological probes. Cluster 


data prefer lower values of these observables deviated at 
more than 2cr level, see e.g. Refs. mm- 

Recently Baryon Acoustic Oscillations (BAO) in the 
Lya forest of BOSS DR11 quasars have been studied at 
redshift z=2.34 nulla- The measured position of the 
BAO peak determines the angular distance, Da(z) and 
expansion rate, H(z). Obtained constraints imply values 
of Da and H that are, respectively, 7% low and 7% high 
compared to the predictions of a flat ACDM cosmological 
model with the best-fit Planck parameters. The signifi¬ 
cance of this discrepancy is approximately 2.5er [12] , 

The tension between CMB based determination of sev¬ 
eral observables by the Planck Collaboration and direct 
low z measurements is intriguing and deserves attention. 
The cause of discrepancy may lie in some calibration er¬ 
rors. On the other hand, it may hint to a deficiency of 
the standard ACDM paradigm. In this paper we show 
that this discrepancy may be resolved if a certain frac¬ 
tion of dark matter is unstable. Decaying Dark Matter 
(DDM) models have been considered previously, see e.g. 
the most recent Refs m H4] , with the stringent con¬ 
straint on DDM decay width T. However, in these papers 
it was assumed that the whole of DM is susceptible to 
the decay, concluding that the decay time must be larger 
than 100 Gyr or so. We instead assume that dark matter 
consists of two fractions, the stable dark matter being 
dominant while a subdominant unstable part decays be¬ 
tween recombination and the present epoch. 

DECAYING DARK MATTER 

Planck constraints. To ensure that our model fits the 
Planck data we accept Planck derived values for all cos¬ 
mological parameters relevant at recombination. In par¬ 
ticular, this means that the sum of initial densities of 
stable and decaying components of dark matter is fixed, 
and after formal redshift to the present moment is deter¬ 
mined by the Planck value io s dm+0Jddm, = 0.1198. In our 
model we vary the initial fraction of decaying component 
in the cosmological mass density 
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FIG. 1. Hubble parameter h as a function of DM decay 
width F for different values of the DDM fraction F. 


FIG. 2. One and two sigma likelihood contours for our model 
parameters. Solid and dashed lines correspond to a dataset 
consisting of JLA sample of SN la and HST measurements of 
h, on top of the best fit Planck model parameters. Addition 
of Planck cluster data results in much narrower shaded area. 


We assume that decay occurs into invisible massless par¬ 
ticles and does not produce too many photons. Alter¬ 
natively, one can consider a scenario when dark mat¬ 
ter consists of two particle species with masses M + p 
and M — p, and the heavier component decays into the 
lighter one with emission of invisible massless particles. 
In this case the dark mass fraction disappearing due to 
decay is equivalent to F = p/M. Throughout the paper 
we normalize the width of the decaying component T to 
km/s/Mpc, i.e. in the same units as Hq. T is another 
independent cosmological parameter in our model which 
we also vary for fitting the data. It is bounded from above 
by the requirement that the unstable fraction does not 
decay substantially before the last scattering to measur¬ 
ably affect the CMB. Hence, we take T < 5000 in which 
range observed CMB spectra are not altered by decays. 

Furthermore, we require that the angular diameter dis¬ 
tance to the last scattering should be the same for all val¬ 
ues of parameters, namely we fix the sound horizon an¬ 
gle 100 *9 S to the Planck value 1.04077. This determines 
Hubble parameter h as a function of F and T and guaran¬ 
tees that derived CMBR spectra in our model are identi¬ 
cal (at high /) to the best fit Plank spectrum for all values 
of parameters. Resulting h as a function of T is shown in 
Fig .[I] for different values of F. Let us remark also that for 
the choice of parameters as in Fig. [l] the age of the Uni¬ 
verse t 0 = |H ( j~ 1 Q)( 1 ' /2 In [(1 + )/(l — U^ 2 )] remains 

nearly the same as predicted by Planck , to « 13.8 Gyr, 
since increasing of Hq is compensated by increasing of 
dark energy fraction Qa- 

Relevant cosmological calculations have been carried 
out using the CLASS Boltzmann code [151 H5] . The pa¬ 
rameter space is explored using the Markov Chain Monte- 
Carlo technique with the Monte Python package [I7j. We 


verified that all CMB spectra are identical at Z > 40. 
At smaller l the spectra somewhat deviate because the 
cosmological constant in our model is typically larger as 
compared to the standard ACDM (we consider spatially 
flat Universe only). However, corresponding changes are 
smaller than the cosmic variance. Therefore we do not 
constrain model parameters using low l Planck data and 
we use supernova data instead. 

Adding supernova and HST constraints. For fitting 
to supernovae observations we use the JLA [18] compi¬ 
lation composed of 740 SN la. This is the largest data 
set to date containing samples from low redshift z « 0.02 
to a large one, z ~ 1.3. The data were obtained from 
the joint analysis of SDSS II and SNLS, improving the 
analysis by means of a recalibration of light curve fitter 
SALT2 and in turn reducing possible systematic errors. 
For ’’standardization” of SN data the linear model for 
the distance modulus p is employed with four nuisance 
parameters in the distance estimates. All necessary data 
for the analysis were retrieved from m- Resulting best 
fit values for all nuisance parameters in our cosmology 
do not differ notably from the values quoted in Ref. T8], 
derived for ACDM. 

We further constrain our model using determination of 
the Hubble parameter with the HST [4j. Resulting one 
and two sigma likelihood contours in the plane of T and 
F are shown in Fig. [2] by solid and dashed lines. We see 
that the base ACDM with T = F = 0 is outside of 2 a 
contours in our model. Derived likelihood for the Hubble 
parameter corresponds h = 0.716 ±0.02 at one a. There¬ 
fore, with a fraction of decaying dark matter the data of 
Planck on CMBR anisotropies, data on supernova, and 
HST data all can be reconciled. 
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FIG. 3. Hubble parameter h(z) (left panel) and angular diameter distance Da (right panel). Model curves are presented for 
fixed T = 2000 and several values of F. Points at non-zero redshift z are the SDSS BAO data. HST measurement at a = 0 is 
also shown with the symbol size comparable to the errorbars. 


DDM and BAO. We now turn to the data on Baryon 
Acoustic Oscillations. The measurement of the charac¬ 
teristic scale of BAO in the correlation function of differ¬ 
ent matter distribution tracers provides a powerful tool 
to probe the cosmic expansion and a convincing method 
for setting cosmological constraints. The BAO peak in 
the correlation function at a redshift z appears at the 
angular separation A 9 = r^/i 1 + z)Da(z), where Da 
is the angular diameter distance and rd = r s (zd ) is the 
sound horizon at the drag redshift, i.e. at the epoch 
when baryons decoupled from photons. BAO feature also 
appears at the redshift separation Az = rd/Du, where 
Djj = c/H{z). Therefore, measurement of the BAO peak 
position at some 2 constrains the combinations of cosmo¬ 
logical parameters that determine Dn/rd and Da/I' d at 
that redshift. 

Recently independent constraints on H rd and Da/ rd 
were obtained using SDSS/BOSS data at z = 0.35 J2H 
EH, z = 0.57 [22l [23] , and z = 2.34 Q2]. These data 
are plotted in Fig. [3] Note that the derived constraints 
for H(z) and Da are not independent, the correlation 
coefficient is 0.5. To avoid cluttering in displaying results 
obtained by different authors at the same redshift, in 
the right panel we plot the results of Refs. HHEIH23] 
for Da/ta, while in the left panel the results of Refs. 
[H E3 Ea for h are presented. Again, the solid line 
corresponds to the ACDM model with the best fit Planck 
measurements. Two other models, F = 0.1 and F = 0.2 
both with T = 2000, are also shown. We see that the data 
systematically deviate from the base ACDM. Thought at 
z < 1 each deviation is about ler (and therefore this is 
not considered as a problem) they all are in the direction 
of the DDM models, except for the result of ref. 23] at 


x = 0.57. 

We repeated procedure of likelihood analysis of previ¬ 
ous subsection with BAO data added. We used BOSS 
BAO likelihoods included into Monte Python pack¬ 
age pi], latest release 2.1. The result is similar to the 
one presented in Fig. [2] but with one and two sigma con¬ 
tours shifted down by about a factor of two. However, 
the result would clearly depend upon the dataset chosen. 
In Ref. E3] DDM model was analyzed with the inclu¬ 
sion of the latest BAO results only, at z = 0.57 [23: and 
x = 2.34 [I2j, with a pessimistic conclusions. In Fig.[3]we 
can see the origin for this conclusion as well. DDM helps 
to ease tension at z = 2.34 both for H(z) and Da, which 
are at the 2.5er level compared to the predictions of the 
base ACDM. However, results of [23] at z = 0.57, which 
are also discrepant at la level, behave differently. While 
DDM is better fit for H{z), it is not so for Da, the latter 
is represented by an upper (red) datapoint at z = 0.57 
in the left panel of Fig. [3] Overall, DDM does not help 
much here. As Ref. 3] describes, at present it is not 
clear whether the discrepancy at z = 2.34 is caused by 
systematics in the Lya BAO measurements (which are 
more complex and less mature than galaxy BAO mea¬ 
surements) or it is an indicator of a new physics. 

DDM and cluster counts. Decaying Dark Matter 
model is capable to resolve tension between the base 
ACDM model and the cluster data as well [25]. This 
is displayed in Fig. [4] in the as and fl m parameter plane. 
The base ACDM corresponds to the error cross marked 
PLANCK CMB. Shaded areas correspond to the pa¬ 
rameter regions allowed (at 2er) by Planck cluster data 
[7, E] and by extended ROSAT-ESO Flux Limited X-ray 
Galaxy Cluster Survey (REFLEX II) [IS]. We should 
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FIG. 4. f lm and ag derived from cluster counts and from 
CMB. Line marked DDM shows trend of these parameters 
when F and F are varied in our model. White circle represents 
a model with F = 0.1 and T = 2000 as an example. 

also note that earlier results obtained in [5;, while in 
agreement with j!0| . are even father away from the 
Planck base ACDM model. In DDM model, when F and 
r are varied, ag and Q m closely follow the line marked 
DDM in Fig.[4]and cross the region allowed by the cluster 
data. White circle on this line represents a model with 
F = 0.1 and T = 2000. With smaller values of F and/or 
r the dot representing a model moves to right, closer to 
the base ACDM model. 

We have added Planck constraints [7] from Sunyaev- 
Zeldovich cluster counts, as (O m /0.27)°' 3 = 0.78 ± 0.01, 
to our likelihood analysis (without BAO). Result is shown 
in Fig. [2] as shaded area, (la area continues actually up 
to F « 0.25 and F « 100 but this is unresolved at the 
scale of this figure.) Now the likelihood of base ACDM 
is vanishingly small compared to a best fit DDM models. 
However, as Planck collaboration concluded on cluster 
counts issue [1]: it is unclear if this tension arise from 
low-level systematics in the astrophysical studies, or rep¬ 
resents the first glimpse of something more important. 
We would say again that the hypothesis of decaying dark 
matter may help to resolve this tension as well. In fact, 
from the joint fit shown in Fig. [2] one can see that the 
issues of Hq and as can be resolved with the same pa¬ 
rameter values of the DDM model. 


CONCLUSIONS 

Cosmological parameters deduced from the Planck 
measurements of the CMB anisotropies with un¬ 
precedented accuracy are at some tension with direct 
astronomical measurements of various parameters at low 
redshifts. We have shown that Planck -inspired ACDM 


cosmology can be reconciled both with HST measure¬ 
ments and with cluster data within the hypotheses of 
Decaying Dark Matter. Joint fit to Planck , supernova, 
HST and Planck cluster data tells that if the dark 
matter decayed between recombination and the present 
time, then the unstable fraction should be about 10 per 
cent at the recombination epoch. Situation with the 
BAO discrepancies is less clear at present and we should 
wait to see in which direction the intrigue will develop. 

Note Added. After our paper has been submit¬ 
ted to arXiv, Ref. j2B] appeared where DDM was 
also suggested as a resolution to the possible tension 
between the CMB and weak lensing determinations of erg. 
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